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Abstract. We present corrections for the change in the apparent scalelengths, central surface brightnesses and axis 
ratios due to the presence of dust in pure disk galaxies, as a function of inclination, central face-on opacity in the 
B-band (r^) and wavelength. The correction factors were derived from simulated images of disk galaxies created 
using geometries for stars and dust which can reproduce the entire spectral energy distribution from the ultraviolet 
(UV) to the Far- infrared (FIR)/submillimeter (submm) and can also account for the observed surface-brightness 
distributions in both the optical/Near-infrared and FIR/submm. We found that dust can significantly affect both 
the scalelength and central surface brightness, inducing variations in the apparent to intrinsic quantities of up to 
50% in scalelength and up to 1.5 magnitudes in central surface brightness. We also identified some astrophysical 
effects for which, although the absolute effect of dust is non-negligible, the predicted variation over a likely range 
in opacity is relatively small, such that an exact knowledge of opacity is not needed. Thus, for a galaxy at a 
typical inclination of 37° and having any > 2, the effect of dust is to increase the scalelength in B relative to 
that in I by a factor of 1.12 ± 0.02 and to change the B-I central colour by 0.36 ± 0.05 magnitudes. Finally we 
use the model to analyse the observed scalelength ratios between B and I for a sample of disk-dominated spiral 
galaxies, finding that the tendency for apparent scalelength to increase with decreasing wavelength is primarily 
due to the effects of dust. 

Key words, galaxies: spiral - galaxies: structure - galaxies: photometry - galaxies: fundamental parameters - ISM 
dust, extinction - radiative transfer 



1. Introduction 

A primary goal of modern studies of star-forming galax- 
ies is to understand how these systems were assembled 
over cosmic time. If the disks of spiral galaxies grow from 
the inside out, as predicted by semi-analytical hierarchical 
models for galaxy formation (e.g. Mo, Mao & White 1998), 
one would predict the stellar populations to be younger 
and have lower metallicity in the outer disk than in the 
inner disk, such that local universe galaxies should be in- 
trinsically larger at the shorter wavelengths where light 
from the young stellar population is more prominent. For 
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the same reason one would expect the intrinsic sizes of spi- 
ral disks to be larger at the current epoch than at higher 
redshift. Observationally, such predictions can be tested in 
two ways. One way is to compare the spatial distribution 
of the constituent stellar populations at different wave- 
lengths, for local universe galaxies. Another way is to look 
for structural differences in galaxies observed at different 
cosmological epochs, at the same rest frame wavelength. 
Both methods require an analysis of the surface-brightness 
distribution of spiral galaxies in the optical and near- 
infrared (NIR) to quantify the distribution of starlight, 
for example by deriving the scalelength of the disk. This 
is done by fitting observed images with models for the 
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surface-brightness distribution of stellar light, whereby 
the disk component is usually specified by an exponential 
distribution. The derived exponential scalclcngths can ci- 
ther be intercompared between different wavelengths for 
local universe galaxies (Elmegreen & Elmegreen 1984, 
Peletier et al. 1994, Evans 1994, De Jong 1996a,b, dc Grijs 
1998, Cunow 1998, 2001, 2004, MacArthur et al. 2003, 
Mollenhoff 2004), or between galaxies at different redshifts 
at a given wavelength (Lilly et al. 1998, Simard et al. 1999, 
Ravindranath et al. 2004, Trujillo & Aguerri 2004, Trujillo 
et al. 2005, Barden et al. 2005). However, the appearance 
of disk galaxies is strongly affected by dust and this effect 
is different at different wavelengths and for different opac- 
ities. This has consequences not only for the derivation of 
the variation of intrinsic scalclcngth with wavelength, but 
also for the variation of intrinsic scalelength with cosmo- 
logical epoch, since the opacity of disk galaxies is expected 
to have been systematically higher in the past ( e.g. Dwek 
1998, Pei et al. 1999). 

The effect of dust on the observed scalelengths and 
central surface brightnesses of disk galaxies has been pre- 
viously modelled by Byun et al. (1994), Evans (1994) 
and Cunow (2001). By means of radiative transfer cal- 
ciilations these works investigated the dependence of 
this effect on star/dust geometry, opacity, inclination 
and wavelength. Recently, a better knowledge of the 
star/dust geometry has been obtained through a joint 
consideration of the direct starlight, emitted in the ul- 
traviolet (UV)/optical/NIR, and of the starlight which 
is re-radiated in the Far-infrared (FIR) /submillimetcr 
(submm). In a series of papers devoted to modelling the 
spectral energy distributions (SEDs) we derived geome- 
tries of the distribution of stellar light and dust that arc 
successful in reproducing not only the observed integrated 
SEDs, but also the observed radial profiles both in the op- 
tical/NIR and FIR/submm (Popcscu ct al. 2000; hereafter 
Paper I, Misiriotis et al. 2001; hereafter Paper II, Popescu 
et al. 2004, see also Popescu & Tuffs 2005 ^). 

In this paper, the fourth in this scries, wc use these 
derived distributions of stars and dust to obtain a new 
quantitative measure of the effect of dust on the ob- 
served photometric parameters in the optical wavebands. 
Furthermore, because it is no longer necessary to explore 
a wide range of star /dust geometries, we are also able to 
systematically explore the full parameter space in opac- 
ity, inclination and wavelength, and tabulate the results 
in a form convenient for the use of the community. We give 
quantitative measures of the change in the apparent scale- 
length, central surface brightnesses and inclination of disk 
galaxies due to the presence of dust.^ All these changes 
are expressed as the ratio of the apparent quantity (i.e. 



^ A simplified version of this geometrical prescription has 
been applied by Misiriotis et al. (2004) to fit the FIR SEDs of 
bright IRAS galaxies. 

^ These corrections are only valid for normal disk galaxies 
and are not applicable to systems with different star/dust ge- 
ometries such as starburst or dwarf galaxies. 



that obtained by fitting images of dusty disks with pure 
exponential disks) to the intrinsic quantity (i.e. that which 
would be obtained in the absence of dust). These correc- 
tions have been derived from a subset of the simulated 
images (those in the optical bands) presented in Tuffs et 
al. (2004; hereafter Paper III). 

In Sect. 2 we give a brief description of the distribu- 
tions of stars and dust used in the simulated images. In 
Sect. 3 wc specify the fitting procedure used to extract ap- 
parent scalelengths, central surface brightnesses and axis 
ratios from the simulated images. These quantities are tab- 
ulated in Sect. 4, where we also describe and explain their 
dependence on opacity, inclination and wavelength due to 
the effect of dust. In Sect. 5 we examine the impact of 
these new results on our ability to derive quantities of 
astrophysical interest from optical observations and give 
in Sect. 6 a specific example of the determination of the 
variation of intrinsic scalelength with wavelength for lo- 
cal universe galaxies. A summary of the paper is given in 
Sect. 7. 

2. The simulated images 

In Paper III we presented simulated images for the diffuse 
component of a spiral galaxy. This diffuse component is 
comprised of a diffuse old stellar population and associated 
dust and a diffuse young stellar population and associated 
dust. Direct evidence for the existence of the diffuse dust 
disk is provided by images of FIR emission in spiral galax- 
ies (for example in M 31 by Haas et al. 1998 and in M 33 
by Hippelein et al. 2003), which clearly show a diffuse disk 
of cold dust emission prominent at 170 /im. FIR measure- 
ments of a statistically unbiased optically selected sample 
of gas-rich galaxies (Tuffs et al. 2002a, b) have shown that 
this diffuse cold dust emission component is ubiquitous 
along the Hubble sequence of late-type galaxies, and car- 
ries the bulk of the dust luminosity for most of these sys- 
tems (Popescu et al. 2002). It is thus responsible for most 
of the attenuation of the stellar light and the associated 
modification of disk brightness and scalelength measured 
in the optical, which is the subject of this paper. 

We know that in reality the apparently diffuse dust 
component may itself have some structure, and that some 
fraction (currently unconstrained by direct observation) 
may be contained in clumps which are externally heated 
and have no embedded sources, the so-called "passive 
clumps" (Popescu & Tuffs 2005). However, providing these 
passive clumps are optically thin, the attenuation charac- 
teristics will be almost identical to those of a completely 
homogeneous distribution (Kuchinski et al. 1998, Pierini 
et al. 2004). If the passive clumps are considered to be op- 
tically thick, then the overall attenuation of a galaxy can 
be reduced by up to 40% (Bianchi et al. 2000a, Misiriotis 
& Bianchi 2002). But, once UV photons can no longer 
penetrate a clump, the photoelectric heating of the gas is 
lost and there is nothing to prevent the clump collapsing 
to form stars. This means that in practice optically thick 
passive clumps should not exist as stable long-lived struc- 
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tures, but rather as transient precursors of star-forming 
clouds. Once the star-formation is underway, the clumps 
will start to contain embedded sources, becoming "active" 
clumps (Popescu & Tuffs 2005). In terms of the global at- 
tenuation of stellar light in galaxies, the main effect of 
active clumps is to locally absorb a fraction of the UV 
output from the embedded massive stars and transform 
it into warm dust emission. This is taken into account by 
the model of Popescu et al. (2000) used in this series of 
papers. Active clumps have such a small filling factor in 
normal spiral galaxies that their effect on the propagation 
of the long-range optical/NIR photons in galaxy disks is 
negligible. Since in this paper we only consider the opti- 
cal/NIR spectral range, it is only necessary here to take 
into account the diffuse dust distribution. 

The diffuse old stellar population has both a "bulge" 
and a "disk" component, whereas the diffuse young stellar 
population resides only in a "thin disk" . Here we only use 
the disk component, since in this paper we are only con- 
cerned with the appearance of the disk in the optical/NIR 
range. The simulations for the disk were performed taking 
the geometry of the diffuse dust to be the superposition 
of the dust in the disk and the dust in the thin disk. In 
other words the attenuated images represent the appear- 
ance of the stellar populations in the disk as seen through 
the dust in the disk and thin disk. The stellar emissivity of 
the disk rj and the extinction coefficients of the dust asso- 
ciated with the disk, k'^^^, and with the thin disk, nl'^f^, 
are described by exponentials: 



'n{X,R,z)='n'^'^''{X,0,0) exp 
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^disk 
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where R and z are the cylindrical coordinates, 77'^'^'^(A, 0, 0) 
is the stellar emissivity at the centre of the disk, /ig'^*^ 



and 



.^disk 



are the scalelength and scaleheight of the stellar 



emissivity of the disk, Kex?(A,0, 0) is the extinction co- 
efficient at the centre of the disk, /i^'^^ and z"^^^^ are the 
scalelength and scaleheight of the dust associated with the 
disk, K*^j**'^(A, 0, 0) is the extinction coefficient at the cen- 
tre of the thin disk and h}^^^^ and z^^^^^ are the scalelength 
and scaleheight of the dust associated with the thin disk. 
The values of these parameters are given in Tables 1 and 
2 from Paper HI. For the old stellar population and as- 
sociated dust the values of these parameters were taken 
from the modelling of Xilouris et al. (1999). It should be 
noted that the thin disk of dust is an approximation to 
the distribution of dust associated with the spiral arms, 
as explained in Paper I, II and III. This approximation 
will be a better representation of spirals having multiple 
spiral arms covering a large fraction of the disk than of 



spirals having a two arm pattern with low coverage of the 
disk or of spirals with a ring structure. This may introduce 
some noise in the attenuation characteristics, in particu- 
lar at the longest NIR wavelength, where the true relative 
geometry of stars and dust may differ the most from that 
assiimed by the model. However it is at these wavelengths 
that attenuation is the least severe. It should also be noted 
that the geometry used by our model is only valid for mas- 
sive (non-dwarf) spiral galaxies, for which Dalcanton et 
al. (2004) have suggested that disk instabilities lead to a 
thin dust layer through collapse and fragmentation of high 
density gas within spiral arms. For low-mass galaxies with 
circular velocities of less than 120km/s Dalcanton et al. 
(2004) have shown that the disk has a different structure, 
with the dust having a larger scaleheight with respect to 
the stars than in more massive, spiral, galaxies. 

The simulated images of the disk were calculated using 
the radiative transfer code of Kylafis & Bahcall (1987), 
which includes anisotropic multiple scattering in which 
the higher orders of scattered light are calculated using 
the method of "scattered intensities" (sec also Kylafis & 
Xilouris 2005). The images have a pixel size (equal to the 
resolution) of 0.0066 of the B-band scalelength /ig'^"^ and 
were sampled every 5 and 10 pixels in the inner and outer 
disk, respectively. They extend out to a radius of 4.63 B- 
band scalelengths /ig"*^, which is equivalent to 3.31 dust 
scalelengths h'^^'^. Examples of simulated disk images are 
shown in Fig. 1 from Paper HI. 

Here we make use of a subset of 280 of the simulated 
disk images presented in Paper III, spanning 7 values of 
total central face-on optical depth in B band Tg, 8 incli- 
nations i and 5 wavelengths. For the sampling in Tg we 
chose the set of vahies 0.1,0.3,0.5,1.0,2.0,4.0,8.0, which 
range from extremely optically thin to moderately opti- 
cally thick cases. For the sampling in inclination we used 
Acos(i) = 0.1, and < (1 — cos(i)) < 0.7, which corre- 
spond to i = 0,26,37,46,53,60,66,73°. The limit of 73° 
in i was imposed since standard methods to derive appar- 
ent exponential scalelength from observed images are no 
longer applicable for higher inclinations. The simulations 
were performed in the standard optical/NIR bands B, V, 
I, J, and K. We also used the corresponding intrinsic im- 
ages of the stellar emissivity (as would be observed in the 
absence of dust). 



3. Derivation of apparent disk scalelengths and 
central surface brightnesses 

3.1. Formalism 

When extracting apparent scalelengths, central surface 
brightnesses and axis ratios from observed (dust attenu- 
ated) images of galaxy disks, the common procedure is 
to fit the observed surface-brightness distribution with 
the brightness distribution corresponding to the projec- 
tion of an inclined thin axisymmetric exponential disk. 
This brightness distribution takes the form: 
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4pp = ^app exp{-^ ) (4) 

where /^pp is the apparent central surface brightness, -Rapp 
is the apparent exponential scalelength, and Qapp is the 
apparent axis ratio b/a oi the inclined disk (equivalent to 
the cosine of the apparent inclination of the disk). The 
coordinate system is orientated such that galaxies have 
their major axis along the x-axis. 

The strategy we adopt in this paper is to use the same 
fitting procedure as used to fit real observations to de- 
rive /appj -Rapp and Qapp from our simulated dust attenu- 
ated disk images. These quantities can then be compared 
with the corresponding intrinsic quantities /q, Rq, and 
Qo, where Iq is the intrinsic face-on surface brightness, 
i?o is the intrinsic exponential scalelength, and Qo is the 
intrinsic inclination of the disk. The intrinsic quantities 
are known as input to our simulations or can be derived 
from the simulated dustless face-on images: Ro is equal to 
the value of the scalelength of the stellar emissivity h'^^^^, 
Qo is cos{i) for a simulation made with inclination i, and 
Iq is the brightness of the central pixel of the simulated 
dustless face-on image. The ratios I'^^^/Iq, Rapp/Ro, and 
Qapp/Qo can then be used to quantify the discrepancy be- 
tween intrinsic and perceived quantities arising from the 
presence of dust. 




Fig. 1. Top: B image of the dusty disk model i = 26°, 
Tg = 8.0 (left) and residuum after subtraction of the pho- 
tometric fit (right). Bottom: B image of the dusty disk 
model i = 73°, Tg — 8.0 (left) and residuum after subtrac- 
tion of the photometric fit (right). For a better visibility an 
alternating intensity table for the disk models was chosen 
which produced the isophote-like images. 



3.2. Technical procedure of the fit 

The simulated disk images containing 700 x 1400 pixels 
and covering half a galaxy were first converted into full 
images of galaxies containing 1399x1400 pixels. The sky- 
value was zero by construction. The maximum of the light 
distribution was determined and was used as the centre of 
the disk fit. 

The two-dimensional surface-brightness function from 
eq. 4 was fitted to the brightness distribution of the simu- 
lated images following Mollenhoff & Heidt (2001). Fitting 
lapp to the two-dimensional surface-brightness distribu- 
tion results in a nonlinear system of equations for the 3 free 
parameters /^pp, i?app and Qapp- A Levenberg-Marquardt 
algorithm was used for the solution of this system (e.g. 
Bevington & Robinson, 1992). 

Several technical steps were necessary as prerequisites 
for the fit: The centre of the galaxy was determined in each 
image and the function from eq. 4 was fitted to the whole 
image. Since there were 280 images to fit, a systematic 
procedure was necessary. We started in K with the face- 
on image and the lowest optical depth {i = 0°, Tg = 0.1). 
A rough guess for the start values of the parameters in the 
initial Levenberg-Marquardt fit was sufficient. The result- 
ing parameters I^pp, Rapp and Qapp were used as start val- 
ues for the next fits {i = 0°, rf = 0.3,0.5,1.0,2.0,4.0,8.0). 
Then the same procedure was repeated with the inclina- 
tions i — 26° until i = 73°. Afterwards we proceeded in 
the same manner for the other filters J, I, V, B. 



Fig. 1 shows two examples of simulated disks and of 
residua after subtraction of our photometric fits to these 
simulated disks. For face-on galaxies the fit is perfect and 
the residuum is practically zero. For inclined galaxies the 
stronger forward scattering induces a relative brighten- 
ing of the near side of the disk, leading to an asymme- 
try around the horizontal axis. Such an asymmetry is not 
considered by our fitting function, therefore the fit to the 
inclined galaxies will be imperfect. 

To show the details of the photometric fits and of the 
residua between the simulated disks and the fits, we plot in 
Fig. 2 cuts parallel and perpendicular to the major axis. 
This is shown for simulations done in the B band, for 
Tg = 1.0 and for three different inclinations, i ~ 0, 46, 73°. 
Due to the central opacity the disk models are somewhat 
rounder in the centre than the corresponding exponen- 
tial fits. This effect is reflected in the central cusp of the 
residuum (solid bottom line) and increases with increasing 
inclination. Nevertheless the exponential fit is a good rep- 
resentation of the simulated brightness distribution. The 
above mentioned asymmetry due to dust scattering is es- 
pecially prominent at high inclination in the cuts perpen- 
dicular to the major axis (bottom right panel of Fig. 2). 

4. Results of the photometric fits 

Tables 1 to 5 show the results of the fitting procedure 
for the B, V, I, J, and K bands. For each table the first 
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Fig. 2. Behaviour of two-dimensional fits and residua for 
simulations done in the B band, for Tg = 1.0, and for 3 
different inclinations: i = (top), i = 46 (centre), and 
i = 73° (bottom). Each panel shows cuts parallel and 
perpendicular to the major axis, taken through the centres 
of the simulated galaxy (solid lines), of the fit (dashed 
lines), and of the residuum (solid lines at bottom). 



two columns contain the input parameters inclination i 
and central face-on central optical depth in the B band 
and the last four columns contain the output parameters 
Qapp, Rapp/Ro, Ilpp/I^ (as defined in Sect. 3.1) and ASS. 
ASB = — 2.51og(F/Fo) is the ratio of the apparent aver- 
age central surface brightness F to the intrinsic average 
central face-on surface brightness Fq, expressed in magni- 
tudes. F is an average taken over an elliptical aperture ^ 
centred on the position of peak brightness in the simulated 
dusty image, with a semi-major axis of i?app/10 and an 
axis ratio of cos(i). i^o is an average taken over a circular 
aperture centred on the position of peak brightness of the 
face-on dustless simulated image, with a radius Rq/1Q. 

^ Wo give the average surface brightness in an aperture (in 
addition to the surface brightness at the centre position of the 
fitted template) because in practice the brightness is blurred 
over a point spread function. Therefore it is useful to include 
a reference aperture which is larger than typical point spread 
functions, but smaller than the scalelength of the disk. 
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Fig. 3. The ratio of the apparent to intrinsic scalelengths 
Rapp/Ro versus Tg, for i = 0, 26, 37, 46, 53, 60, 66, 73°, and 
for the wavebands B, V, I, J and K. 

Here we should mention that the results given here 
are, as described in Sect. 3.2, obtained by fixing the co- 
ordinates of the centre of the fit to the position of max- 
imum observed brightness. Other fitting procedures such 
as GIM2D (Simard 1998, Marleau & Simard 1999) or 
GALFIT (Peng et al. 2002) constrain the coordinates of 
the centre from the moments of the two-dimensional in- 
tensity distribution, which result in a small shift along the 
minor axis of the fitted centre for opaque high-inclination 
disks. The effect of this shift on the derivation of the cen- 
tral brightness is at most 3%. The effect on the derivation 
of the scalelength is negligible, since it depends only to 
second order on the shift. 

In the following subsections we describe the general 
behaviour of the structural parameters listed in Tables 1- 
5 as a function of , i and wavelength. 

4.1. Change of the exponential scalelengths 

In Fig. 3 we show the dependence of Re,pp/Ro on Tg, for 
different inclinations and wavelengths. Because disks are 
always more attenuated in their central regions than in 
their outer regions, R^pp/Ro is always greater than unity 
and increases monotonically with Tg. In some cases this 
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increase is linear, whereas in others the increase flattens 
off towards higher vahics of Tg . The linear behaviour hap- 
pens when the disk is optically thin along the lines of sight 
towards most of the projected radii. As Tg increases, an 
increasing area of the inner disk becomes optically thick, 
producing a slowing-down in the increase with Tg, until 
an asymptote is approached when the disk becomes opti- 
cally thick along the lines of sight over the whole extent 
of the disk. At longer wavelengths, e.g. in the K band, 
such an asymptote is never reached and the behaviour of 
the curves remains in the linear regime, since even at the 
highest Tg and inclinations considered, the disk remains 
optically thin along most of the lines of sight. By contrast, 
at the shorter wavelengths, for instance in the B band, the 
asymptote is approached for high Tg and inclinations, as 
the disk becomes optically thick along most of the lines of 
sight. Furthermore, at these wavelengths the curves enter 
the non-linear regime at all inclinations, since some part 
of the disk is always optically thick along the line of sight. 
At intermediate wavelengths, such as in the I band, the 
variation with Tg remains in the linear regime at lower 
inclinations, and only at high inclination is a non-linear 
behaviour seen, though without getting near the asymp- 
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Fig. 5. The ratio of the apparent to intrin- 
sic scalelengths Rapp/Ro versus waveband, for 
i = 0,26,37,46,53,60,66,73°, and for t| = 
0.1,0.3,0.5,1.0,2.0,4.0,8.0 



tote. Another consequence of the existence of the asymp- 
tote is that the curves for different inclinations tend to 
bunch together at high Tg and shorter wavelengths, as 
the asymptote is approached. Thus, in B, V and I bands, 
the highest spread between the curves at different incli- 
nations is reached at intermediate Tg, since, by necessity, 
the curves must also be bunched together at low Tg . In K 
band the curves never enter the non-linear regime so the 
maximum spread happens at the highest Tg. 

In Fig. 4 we plot the same information as in Fig. 3, but 
as the variation of Rapp/Ro with inclination, for different 
Tg and wavelengths. The curves appear flat when the disk 
is either predominantly optically thin at all inclinations 
(for example in K band at t| < 2, where Ra,pp/Ro is close 
to unity) or when the disk is predominantly optically thick 
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Fig. 6. The change in the ratio of the scalelengths in 
B and I bands due to the effects of dust, plotted as 



^app(-g) 
-Kapp(-f) 



versus Tg. Different curves are plot- 



ted for a range of inclinations from z = to i = 73. The 
curves for intermediate inclinations {i = 46, i = 53 and 
i = 60) are plotted with dotted lines for clarity. 
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Fig. 7. Comparison between the apparent ratios of scale- 
length in B and I predicted by our work (solid lines) and 
by the work of Cunow et al. (2001) (dashed lines). The 
curves are for Tg = 0.1, 1, 4 and 8 (interpolated in the 
case of the Cunow curves). The pairwise corresponding 
curves are marked by 4 different symbols. 



at all inclinations (for example in B band at Tg =8, where 
Rapp/Ro approaches the asymptote of 1.5). 

In Fig. 5 we plot Rapp/Ro versus wavelength, for dif- 
ferent Tg and inclinations. The variation with wavelength 
also shows a transition between a linear behaviour in an 
optically thin regime (for t| < 0.5) and a non- linear ap- 
proach to an asymptote at Tg = 8. 

In order to see what is the relative change in the scale- 
length between different wavelengths just due to the ef- 
fect of dust, one can look at the variation of the quan- 



tity ^^ppi^^) X -^0(^2) 

-fi'app(A2) i?o(Al) 



with Tg . As an example, we plot 



in Fig. 6 this quantity for the wavelengths B and I. At 
all inclinations, the curves show a similar basic shape. 
The curves rise steeply at low Tg, reach a maximum and 
then decline slowly. This behaviour is compressed over a 
smaller range of Tg for higher inclinations. Because of this, 

-^app (B) ^ ^o(£) increases with inclination at low t£, 

-rtapp(,-/j rCo[±!) 

but decreases with inclination at high Tg, with the curves 

f f 
crossing at intermediate Tg. In other words, for a low Tg, 

in order to convert the ratio between the apparent scale- 
length in B and I to the corresponding intrinsic ratio, one 
needs to apply a correction factor which is larger for high 
inclinations than for low inclinations, whereas at high Tg 
the opposite is true. An interesting feature of these curves 
is that they are relatively flat for Tg > 2. 

Simulations of the observed ratios of scalelength in B 
and I as a function of inclination and Tg have also been 



performed by Cunow et al. (2001). The closest compari- 
son possible is that with model (k) of Cunow, which differs 
from our model only through the absence of the dust disk 
associated with the young stellar population. In Fig. 7 we 
have plotted the ratios of apparent scalelength in B and I 
as a function of apparent cUipticity {1 — b/a) for both our 
model and the model from Fig. 7 of Cunow et al. (2001), 
where the ratios from Cunow have been interpolated to 
match the opacities given in our tables. Overall the curves 
of Cunow behave as if their model is more optically thin for 
the same Tg, which is the opposite behaviour one would 
expect from the differences in the geometrical distribu- 
tions between the models. However the procedure to derive 
scalelength followed by Cunow et al. only uses data from 
the outer parts of the disk having a linear profile, whereas 
our method uses data from the whole disk. This obviously 
biases the results from Cunow towards more optically thin 
solutions and this effect dominates over the differences in- 
duced by the different geometrical distributions. 

Another work on predicting apparent scalelength ra- 
tios between B and I is that of Byun et al. (1994). These 
results are very similar to those of Cunow et al. (2001) (see 
Fig. 8 of Cunow et al.), which may reflect the fact that 
Byun et al. also use extrapolation from the outer disk in 
deriving photometric parameters. 

4.2. Change of the central surface brightness 

In Fig. 8 we show the dependence of ASB on Tg, for differ- 
ent inclinations and wavelengths. An interesting feature of 
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Fig. 8. The ratio of the apparent average central surface 
brightness to the intrinsic average central face-on surface 

brightness expressed in magnitudes, ASB, versus Tg, for 
i = 0,26,37,46,53,60,66,73°, and for the wavebands B, 
V, I, J, and K. 
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Fig. 9. The ratio of the apparent average central sur- 
face brightness to the intrinsic average central face-on 
surface brightness expressed in magnitudes, ASB, versus 
l-cos{i), iorrl = 0.1, 0.3, 0.5, 1.0, 2.0, 4.0, 8.0 and for the 
wavebands B, V, I, J, and K. 



these ciirvcs is that they pass through both positive and 
negative values of ASB: the centre of the disk can appear 
either dimmer or brighter than the face-on dustless disk. 
This is because the projected surface density of stars visi- 
ble along the line of sight can be either greater or smaller 
than the surface density of stars seen towards the centre 
of a face-on dustless disk. The projected surface density 
of stars can be increased by viewing the disk at higher in- 
clinations and can be decreased by attenuating the stellar 
light with dust. In the K band the central part of the disk 
appears brighter than the centre of the face-on dustless 
disk at almost all inclinations and Tg, because the line of 
sight through the centre is almost always transparent. On 
the contrary, in the B band the opposite is true, since the 
line of sight through the centre is almost always opaque. 

The form of the curves in Fig. 8 can be explained 
in terms of the two competing factors viewing angle and 
opacity: the brightening of the ^B^pp with increasing in- 
clination and the dimming of the SB^^^ with increasing 

Tg. The former effect is independent of wavelength, since 
the line of sight of stars is always the same for a certain 



inclination^ This can be seen from the wavelength inde- 
pendence of the spread in ASB from to -1 magnitude 
for Tg =0. The latter effect (the dimming due to opacity) 
is strongly dependent on wavelength, since, for a certain 
Tg, the depth into which stars can be seen along the line 
of sight decreases with decreasing wavelength. The dim- 
ming of SB^pp with increasing Tg is almost negligible in 
the K band at inclinations close to face-on, when the line 
of sight through the centre is almost transparent. On the 
contrary, in the B band this effect becomes strong at all 
inclinations. 

Another feature of these curves is the transition be- 
tween a linear dependence of ASB on Tg when the lines of 
sight through the centre arc optically thin (e.g. in the K 
band at all inclinations), to a non-linear regime for opti- 
cally thick lines of sight (e.g. in the B band at all inclina- 
tions). At intermediate wavelengths (e.g. in the J band) 

^ Formally there is a very small dependence on wavelength 
due to the fact that the scalelength of the stars varies with 
wavelength in our model. In practice, however, this effect is 
negligible over the range of inclinations considered here and 
cannot be seen in the plot. 
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Fig. 10. The ratio of tlie apparent average central sur- 
face brightness to the intrinsic average central face-on 
surface brightness expressed in magnitudes, ASB, ver- 
sus waveband, for i = 0, 26, 37, 46, 53, 60, 66, 73°, and for 
t| = 0.1, 0.3, 0.5, 1.0, 2.0, 4.0, 8.0. 



one sees a transition between these two regimes: the upper 
curves (close to face-on) are in the linear regime, whereas 
the lower curves (high inclinations) are non-linear. 

For a given (non-zero) Tg, the dispersion between the 
curves in Fig. 8 for low and high inclinations varies sys- 
tematically with wavelength. For the case of Tg = 4 one 
can see that the dispersion has a maximum in the K band, 
then it decreases towards I band, after which it increases 
again towards B band. At first glance this behaviour is not 
very intuitive. If the stars and dust had an identical dis- 
tribution, then one would expect a monotonic decrease in 
dispersion when moving from transparent central disks to 
opaque central disks, until the lines of sight would become 
completely opaque at all inclinations and the dispersion 
would tend to zero (since the brightness of an optically 
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Fig. 11. The change in the B-I central colour due to the ef- 
fects of dust, plotted as ASB{B) — AS B{I)+ const versus 
Tg . Different curves are plotted for a range of inclinations 
from I = to z = 73. The curves for intermediate inclina- 
tions (i = 46, i = 53 and i = 60) are plotted with dotted 
lines for clarity. 



thick body is independent of the viewing angle). The fact 
that the behaviour is more complex than this is a conse- 
quence of the non-identical distribution of stars and dust 
in the vertical direction. In our model the stars and dust 
are distributed in exponential disks, with the stars hav- 
ing scalehcights larger than those of the dust. As already 
remarked by Disney et al. (1989), "even in a highly ob- 
scured galaxy, the upper layers of unobscured stars will 
apparently increase its mean surface brightness as it is in- 
clined, mimicking the behaviour of an optically thin sys- 
tem". Thus, at the extreme of very high optical depth 
one approaches the same dispersion as for the very low 
optical depth, since the variation in the surface bright- 
ness between low and high inclination is only produced 
by the stars above the (exponential) dust layer. At inter- 
mediate optical depths the dispersion over inclination is 
lower, because a larger fraction of this dispersion comes 
from stars well mixed with the dust which have a smaller 
variation in the surface brightness between low and high 
inclination. In our model, there will always be an increase 
in the surface brightness with increasing inclination, even 
for optically thick disks. This can be seen even better in 
Fig. 9, where we plot the same information as in Fig. 8, 
but as the variation of ASB with inclination, for different 
Tg and wavelengths. We also plotted the variation ASB 
with wavelength (Fig. 10), for different Tg and inclina- 
tions, to better capture the transition between the almost 



independence of ASB on wavelength at low Tg to a strong 
dependence on wavelength at high Tg . 
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In order to see what is the change in the central colour 
of the disk due to dust, one can look at the variation of 
ASB{Xi) - ASB{X2) with r^. As for the analogous plot 
for the variation of scalelength, we plot in Fig. 11 this 
quantity for the wavebands B and I, ASB{B) — ASB{I), 
which is equivalent to a plot of SB(B) - SB (I) shifted by a 
constant value, such that the curves go through the origin. 
These curves, therefore, give the colour correction that 
needs to be applied to correct apparent colour to intrinsic 
colour at the centre of the disk. The curves show the same 
qualitative behaviour as those for the relative change in 
scalelength, namely that they rise steeply at low Tg, reach 
a maximum and then decline slowly. We note again the 
relatively flat behaviour of the curves for > 2. 



4.3. Change of the axis ratios 

In Fig. 12 wc show the dependence of the ratio of the 
apparent to the intrinsic axis ratio, Qapp/Qo on r^, for 
different inclinations and wavelengths. One can see that 
this ratio can be cither greater or less than one, which 
means that the galaxy can appear less or more inclined 
than it is in reality. In other words the galaxy can ap- 
pear rounder or flatter than the shape of an infinitely 
thin dustless disk at the true inclination of the galaxy. 
This behaviour is attributed to two eff'ects that act in the 
opposite sense. Firstly, the vertical distribution of stars 
makes the galaxy appear rounder than an infinitely thin 
dustless disk. Secondly, the presence of dust makes the 
galaxy appear flatter than an infinitely thin diistlcss disk. 
The first effect is obvious, but the second needs some ex- 
plaining. As one increases the line of sight optical depth, 
cither by increasing the inclination or the Tg , the observed 
light becomes dominated by photons originating from an 
increasingly thin stellar layer on the side of the dust disk 
nearest to the observer, due to absorption by dust of pho- 
tons originating from below this layer. 

The eff'ect of the vertical distribution of stars is best 
seen by looking at the inclination dependence of Qapp/Qo 
for Tg close to zero. For example in the K band, Qapp/ Qo 
increases monotonically with inclination up to ~ 1.13 for 
i = 73°. 

The effect of opacity is more complex. One feature is 
the decrease of Qapp/Qo with increasing t^, for all in- 
clinations and wavelengths. This effect is weak in the K 
band where opacity is small, and is strong in the B band 
where opacity is large. A further aspect of the effect of 
dust is that at high opacity, and except for inclinations 
near to edge-on, the dependence of Qapp/Qo on inclina- 
tion is the opposite of that predicted by the model at low 
opacity. Thus, in B-band at Tg =8, the general trend is 
for Qapp/Qo to decrease with increasing inclination^ over 
the range i = to i = 66, as one progressively sees a 
thinner and thinner layer of stars above the dust layer, 
making the galaxy appear flatter with increasing inclina- 
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Fig. 12. The ratio of the apparent to intrinsic axis ratio 



o/Qo 



versus t, 



B' 



for 



0, 37, 53, 66, 73° and for the 



* within the numerical noise induced by the finite numbers 
of directions for which scattered light is calculated. 



wavebands B, V, I, J, and K. 



tion. Suddenly, at i = 73, this trend is reversed, as one 
starts to see stars both above and below the dust layer, 
and a big step is made to a rounder shape. The differ- 
ent behaviour of Qapp/Qo on inclination at low and high 
Tg = 8 means that the curves cross at intermediate Tg. 

5. Discussion 

Having established the dependence on inclination, face- 
on opacity and wavelength of the corrections needed to 
convert the apparent photometric parameters scalelength, 
central surface brightness and axis ratios into the corre- 
sponding intrinsic quantities, we can now discuss the im- 
plications of these findings. But before doing this, we cau- 
tion that the corrections given in this paper are valid only 
for disk-dominated galaxies. In reality, many galaxies have 
prominent bulges and therefore the apparent photometric 
parameters of the disk (including scalelength, central sur- 
face brightness and total magnitudes) have to be found in 
a simultaneous fit of geometrical functions for the bulge 
and the disk. However, the apparent photometric param- 
eters derived from this bulge-to-disk decomposition will 
be a function of both opacity and bidgc-to-disk ratio. In 
other words, the apparent scalelength and central surface 
brightness that one derives for a disk galaxy with a given 
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opacity and inclination in the presence of a bulge will in 
general be different from the apparent scalelength that 
one derives for the same disk (with the same opacity and 
inclination) without a bulge. Therefore one should restrict 
the analysis to galaxies having small bulges, where such 
effects can be ignored. 

One should also caution that the corrections for dust 
given in this paper are calculated assuming a small wave- 
length dependence of the ratio of the intrinsic scalelength 
of the stars to that of the dust, as specified in Tuffs et 
al. (2004). Thus, strictly speaking, the tabulated correc- 
tions arc valid only for this geometry. In practice, however, 
these corrections can be used also for systems having a dif- 
ferent variation in the ratio of the intrinsic scalelengths. 
Comparison between Fig. 7b and h from Cunow ct al. 
(2001) show that, if one were to make an error of 25% in 
the assumed intrinsic scalelength of stars relative to the 
dust, the corresponding error in the effect of dust would 
be at most a few percent for typical opacities and incli- 
nations less than 70°. Only for > 7 and inclination 
greater than 70° does the error exceed 10%. 

After these words of caution we can now examine 
the impact of these new results on our ability to de- 
rive quantities of astrophysical interest from optical ob- 
servations of disk- dominated galaxies. In principle, if one 
knew the face-on opacity, one could use the correction 
factors to analyse the intrinsic properties of galaxies on 
an object to object basis. In practice, however, the face- 
on opacity of an individual galaxy can only be reliably 
determined by sclf-consistently modelling its combined 
UV/optical/FIR/submm output, which is at present only 
rarely possible, since the vast majority of galaxies has not 
yet been measured in the FIR/submm regime. For large 
statistical samples, though, one can derive a characteristic 
value for the face-on opacity of a population of systems 
from the optical measurements themselves by analysing 
the inclination dependence of the photometric parameters 
described here. 

One can also use our tabulated corrections for the in- 
vestigation of astrophysical effects for which, although the 
absolute effect of dust is non-negligible, the predicted vari- 
ation over a likely range in opacity is relatively small, such 
that an exact knowledge of opacity is not needed. An ex- 
ample is the variation of the scalelength between B and I 
due to dust, which is small, provided that Tg > 2 (galaxies 
are not too optically thin).^ Indeed, inspection of Fig. 6 
shows that for a galaxy at a typical inclination of 37° 
and having any Tg > 2, the effect of dust is to increase 
the scalelength in B relative to that in I by a factor of 
1.12 ± 0.02. The scatter in this ratio is so small that one 
can compare the intrinsic scalelength of the stellar popu- 
lations in B and I by simply applying a single correction 
factor for dust for any spiral galaxy at this inclination. 
Similarly, other correction factors can be derived from our 

® In our detailed SED modelling of individual galaxies given 
in Paper I and II we showed that this range in opacity is quite 
likely to characterise spiral galaxies. 
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tables for other inclinations and other waveband combi- 
nations. 



Correspondingly, one can also compare the surface 
brightnesses between B and I to get information about 
the intrinsic colour of the central and most dust affected 
regions of the disks, even without detailed knowledge of 
Tg. Inspection of Fig. 11 shows that for a galaxy at an 
inclination of 37° and for Tg > 2, the change in B-I cen- 
tral colour due to dust is 0.36 ± 0.05 magnitudes. Again, 
the scatter in the central colour is small enough that one 
can look at the intrinsic central colour between B and I 
by simply applying a single correction factor for dust for 
any spiral galaxy at this inclination. 



Lack of knowledge of opacity may be somewhat more 
problematic when one attempts to compare the proper- 
ties of dusty galaxies in the local and distant universe at 
a fixed restwavelength. When studying the evolution of 
disk sizes with redshift (e.g. by comparing the scalelength 
of low {z ~ 0) and high (z ~ 1) redshifted galaxies at a 
common restwavelength) one must correct for the increase 
in the apparent disk scalelength with increasing redshift, 
stemming from the expected increase in opacity. If, for 
example, the opacity increased from 4 at 2 = to 8 or 
more at z = 1, the apparent scalelength in the B band 
of a galaxy seen at i = 37° would increase by a factor 
of at least 1.10 due to dust. This factor would rise to at 
least 1.24 for an increase in opacity from 2 to 8 or more. 
The increase in the apparent scalelength with increasing 
redshift due to the effect of dust goes in the opposite di- 
rection to the predicted evolution in intrinsic disk size, 
and in general would need to be taken into account in 
quantitative analyses. Ultimately one would need to refine 
the technique to take into account the possible evolution 
of the dust/star geometry'' when making corrections for 
high-redshifted galaxies. This will become possible when 
spatially resolved dust emission data from the rest-frame 
far-infrared is available for such systems, thereby permit- 
ting an analysis of the star/dust geometry to be performed 
analogous to that performed for local universe galaxies. 



Our tabulated corrections for dust will also allow the 
analysis of the correlation between the central face-on 
surface brightness versus scalelength, as investigated by 
Graham (2001) in the K band where the galaxies are opti- 
cally thin, to be extended to the shorter wavebands where 
the galaxies may be optically thick. As a by-product of 
such an investigation, we note that one could obtain a 
statistical determination of Tg from the inclination depen- 
dence of both central surface-brightness and scalelength. 
As previously cautioned, however, this analysis has for the 
moment to be restricted to disk-dominated galaxies. 



12 



C. MoUenhoff et al.: Modelling the spectral energy distribution of galaxies. IV 



1.6 



1.4 - 



a 
a 

IS 



a 

(0 



1.2 



0.8 



0.6 





o 
o 


o 

o 






° Oo°° 




o o 


- - -o- - _ 




o 

8 





1 


o 



_0 T=2 
T=4 

- T=8 

- T=0.1- 



_L 



0.2 



0.4 

- b/a 



0.6 



0.8 



Fig. 13. Observed scalelength ratios between B and I, 
R (B) 

appv ; ^ versus observed axis ratios 1 — b/a. The lines 

represent the model predictions for the same relation due 
only to the effect of dust, which is equivalent to the predic- 
tion of apparent scalelength ratios for dusty disks having 
the same intrinsic scalelength in B and I. The lines cor- 
respond to Tg = 0.1,2,4,8. A typical error bar of 7% is 
plotted in the lower right hand corner. 



6. Application to local universe disk galaxies 

As mentioned in the introduction, one of the motivations 
of this work is to quantify the change in the intrinsic scale- 
length of different stellar populations in local universe disk 
galaxies. Here we give an example of the application of 
our model to such an investigation. For this we used pub- 
lished apparent disk scalclcngths in B and I, restricting 
the data to galaxies which have measured bulge-to-disk 
ratios of less than 0.2. There are two reasons for this selec- 
tion criterion. Firstly, as mentioned in the discussion, one 
needs to restrict the analysis to disks with small bulges, 
since the apparent scalelength derived for dusty disks with 
large bulges will diverge from the apparent scalelength of 
pure dusty disks. Secondly, for systems with small bulges 
MacArthur et al. (2003) showed that results from one- 
dimensional photometry are equivalent to those derived 
from two-dimensional fits (as required by this work). The 
resulting sample is comprised of 61 galaxies taken from 
MoUenhoff (2004), MacArthur et al. (2003), Prieto et al. 
(2001) and Pompei & Natah (1997). Since MacArthur et 
al. had measurements in B V R H, their R band scale- 
lengths were converted to I band scalclcngths by multi- 
plying them by a factor of 0.981. This factor was obtained 



^ For example, depending on the assumed star-formation his- 
tory, it may be necessary to invoke different star/dust geome- 
tries to explain the observed colours of distant disk galaxies 
(Pierini et al. 2005). 



by linearly interpolating between the average scalelength 
in R and H, as given in Table 3 of MacArthur et al. . 
Despite our selection of disk-dominated galaxies and the 
known tendency of such systems to preferentially popu- 
late the intermediate- and low-mass end of the mass dis- 
tribution of nearby, late-type galaxies (e.g. Pierini & Tuffs 
1999), the galaxies used in this sample are all bright spiral 
galaxies with absolute magnitudes My = -19.2 to -22.3. 
Thus they have geometries of stars and dust appropriate 
to giant spirals, which are well described by our model. 

In Fig. 13 we plot observed scalelength ratios between 

R (B) 

B and I, fT\ : versus axis ratios 1 — b/a. All but three 

of the 61 data points lie above unity, which means that 
almost all galaxies have apparent scalclcngths which are 
larger in B than in I. The data also show a scatter, which 
is consistent with the random errors estimated for the 
measurements. * The mean and the standard deviation in 
i?app(-B)/i?app(-f) is 1.156 and 0.094, respectively. In order 
to see whether these observed values can be accounted for 
by the effect of dust alone, we compare them with model 
predictions for Tg > 2 (see Fig. 13) which is a reason- 
able range in to expect for disk galaxies. For each of 
the model predictions we calculated the corresponding ra- 
tio of intrinsic scalclcngths at B and 1 which best fit the 
data by shifting the lines vertically in the figure until an 
equal number of data points appeared above and below 
the line. This process yielded average ratios of intrinsic 
scalclcngths of 1.057, 1.019 and 1.037, for t| = 2,4,8, 
respectively, and an overall range (±1(T for the sample) 
of intrinsic scalelength ratios from 0.93 to 1.15. Because 
the deviation from unity of the averaged ratios is so much 
less than the spread in the individual data points, we con- 
clude that, on average, and within the available statis- 
tics, there is almost no systematic increase in the intrin- 
sic scalelength with decreasing wavelength and that the 
observed predominance of galaxies with ratios of appar- 
ent scalelengths greater than unity is mainly due to the 
effects of dust. The data would only admit a small vari- 
ation of the intrinsic stellar scalelength between B and I 
(less than 1.15) if < 2. However combined analysis of 
the UV/optical and FIR/submm have generally yielded 
Tg < 2 (Bianchi et al. 2000b, Paper I, II, Misiriotis et al. 
2004, Meijerink et al. 2005). 

Much bigger statistical samples will be needed to mea- 
sure systematic trends in intrinsic scalelengths with wave- 
length, if they exist. Finally we note that the scatter in 
the observed ratio of scalelengths between B and I (with 
a = 0.094) is comparable to the estimated errors in this 
quantity for individual galaxies and can therefore primar- 
ily be ascribed to measurement errors. 



* The errors of a single scalelength determination were esti- 
mated to be 5% in MoUenhoff, ~ 3% in Cunow and ^ 5% 
in MacArthur et al.. Thus the error of the quotients in Fig. 13 
can be estimated to be ~ 7%. A corresponding error bar is 
plotted there. 
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Our conclusion that dust is mainly responsible for the 
larger scalelength observed in the B band compared to 
I band can be contrasted with the results from Cunow 
(2001, 2004), who deduced that dust alone was not suffi- 
cient to explain this effect. Here we note that the different 
conclusions can primarily be attributed to differences in 
the data rather than differences in the model. In fact, as 
noted by Cunow (2004), her data for normal galaxies, both 
in terms of the absolute values of the scalelength ratios 
and the variation of these ratios with inclination, cannot 
be fitted by any model. As shown in Fig. 13, our data is 
reasonably well fitted by our model in both these respects. 

Another study in the literature that has addressed the 
wavelength dependence of scalelength (this time between 
B and K bands) is that of Peleticr ct al. (1994). Here again, 
on average, the scalelength ratios have much larger values 
than in our sample, even taken into account the differ- 
ent wavelength used. Although Pelctier ct al. interpreted 
their results as being due to the effect of dust on the ap- 
parent scalelengths, such large effects are not predicted by 
any dust model. Bearing in mind the inclusion of galaxies 
with bulges in the sample of Peletier et al., such effects 
are, in our view, more likely to be produced by the effect 
of dust on the bulge-disk decomposition itself. This is be- 
cause the use of simple templates for dustless bulge and 
disk in bulge-disk decomposition would lead to systematic 
bias in the derived photometric parameters, including ap- 
parent disk scalelength. 

7. Conclusions and summary 

We have fitted the simulated images of dusty disk galaxies 
presented in Paper III with a template brightness distri- 
bution corresponding to an inclined infinitely thin dust- 
less exponential disk to obtain apparent scalelengths, cen- 
tral surface-brightness distributions and axis ratios. These 
are the apparent photometric quantities that an observer 
would extract from observed images of galaxies. Using the 
prior knowledge of the corresponding intrinsic quantities 
input in the simulations, we were able to derive the cor- 
rection factors (listed in Tables 1-5 of this paper) for the 
conversion of the apparent to intrinsic quantities due to 
dust, as a function of inclination i, central face-on optical 
depth in B band Tg, and wavelength. 

The apparent to intrinsic scalelength ratio is always 
greater than unity and can vary up to 50%. The appar- 
ent to intrinsic central surface-brightness ratio expressed 
in magnitudes can be either positive or negative, depend- 
ing on whether the dimming due to dust is more or less 
important than the brightening due to the increase in the 
column density of stars induced by inclining the disk. This 
ratio can change up to 1.5 magnitudes due to the effect 
of dust. The ratio of the apparent to intrinsic axis ratio is 
not strongly affected by dust, but has an opposite depen- 
dence on inclination according to whether the lines of sight 
through the disk are predominantly optically thin or opti- 
cally thick. In the former case the vertical distribution of 
stars makes the galaxy appear progressively rounder with 



increasing inclination than an infinitely thin disk, whereas 
in the latter case the opposite trend is seen, because an 
increasing proportion of the observed light originates from 
a thin layer of stars above the dust. 

Assuming that the basic geometry of dust and stars in 
local universe spiral galaxies also applies for higher red- 
shift spiral galaxies, our tabulated corrections can be used 
to correct for the increase in the apparent disk scalelength 
with increasing redshift due to the expected increase in 
opacity. This will allow the intrinsic evolution of disk sizes 
with cosmological epoch to be investigated. As an exam- 
ple, we show that for a possible variation in opacity be- 
tween 2 at z = to 8 or more aX z = 1, the apparent 
scalelength in the B band of a galaxy seen at i = 37° 
would increase by a factor of at least 1.24 due to dust. 

We used our model to analyse the distribution of ob- 
served scalelength ratios between B and I for a sample 
of disk-dominated spiral galaxies. We found that the pre- 
dominance of galaxies with larger apparent scalelength in 
B than in I is primarily due to the effects of dust. 
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Table 1. Structural parameters of the disks in B. 
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Table 4. Structural parameters of the disks in J. Table 5. Structural parameters of the 
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